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SUMMARY 

The presence of mesoscopic heterogeneities in fluid-saturated porous rocks can produce 
measurable seismoelectric signals due to wave-induced fluid flow between regions of dif¬ 
fering compressibility. The dependence of these signals on the petrophysical and struc¬ 
tural characteristics of the probed rock mass remains largely unexplored. In this work, 
we derive an analytical solution to describe the seismoelectric response of a rock sample, 
containing a horizontal layer at its center, that is subjected to an oscillatory compress¬ 
ibility test. We then adapt this general solution to compute the seismoelectric signature 
of a particular case related to a sample that is permeated by a horizontal fracture lo¬ 
cated at its center. Analyses of the general and particular solutions are performed to study 
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the impact of different petrophysical and structural parameters on the seismoelectric re¬ 
sponse. We find that the amplitude of the seismoelectric signal is directly proportional 
to the applied stress, to the Skempton coefficient contrast between the host rock and the 
layer, and to a weighted average of the effective excess charge of the two materials. Our 
results also demonstrate that the frequency at which the maximum electrical potential 
amplitude prevails does not depend on the applied stress or the Skempton coefficient con¬ 
trast. In presence of strong permeability variations, this frequency is rather controlled by 
the permeability and thickness of the less permeable material. The results of this study 
thus indicate that seismoelectric measurements can potentially be used to estimate key 
mechanical and hydraulic rock properties of mesoscopic heterogeneities, such as com¬ 
pressibility, permeability, and fracture compliance. 

Key words: Hydrogeophysics - Electrical properties - Fracture and flow - Permeability 
and porosity. 


1 INTRODUCTION 


Free electrical charges at the mineral surfaces of fluid-saturated porous rocks are responsible for an 
electrical double layer (EDL) in the pore space surrounding the solid grains. The EDL is characterized 
by an electrical excess charge that counter-balances the surface charges (Hunter 19811. Pore fluid flow 
in response to a propagating seismic wave (Biot 1956a|b ) exerts a drag on this excess charge, thus 
resulting in an electrical source current generally referred to as streaming current. This in turn results 
in an electrical potential distribution, which can be measured in laboratory and field experiments. 
The thus generated electric field depends on a range of important petrophysical parameters, such as 
the porosity, the permeability, and the type of saturating fluid. This has fostered an increased interest 


in this physical process, commonly referred to as seismoelectric conversion (e.g., Thompson & Gist 
T993i rJouniaux [20TTj ). 


Pride (1994) developed the theoretical basis of the seismoelectric phenomenon by coupling Biot’s 


(1956a, 1956b) and Maxwell's equations through a volume-averaging approach. This theory predicts 
two kinds of seismoelectric conversions: (1) the coseismic field and (2) the interface response. The 
coseismic field is a consequence of the wavelength-scale relative fluid flow associated with the passing 
seismic wavefield, which in turn generates a streaming current and thus an electric field. This field 
travels with the seismic wave, even in the case of homogeneous media, but is largely limited to the 
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wave’s support volume (Pride & Haartsen 1996). Conversely, the interface response occurs when a 
seismic perturbation strikes on an interface in terms of mechanical or electrical properties. In this case, 
a variation in the streaming current distribution arises, breaking the symmetry in the charge separation. 
This generates an electromagnetic (EM) signal that propagates independently of the seismic wave and 
can be measured outside the support volume of the wave, thus allowing for the remote detection of 
geological interfaces. These signals are highly sensitive to the fluid pressure gradients in the vicinity of 
the interface. Thus, proper modeling of wave conversions at interfaces, and in particular of Biot slow 


waves, is critical for the accurate modeling of seismolectric interface responses (Pride & Garambois 


2002 ). 


The propagation of seismic waves through a medium containing mesoscopic heterogeneities, that 
is, heterogeneities having sizes larger than the typical pore scale but smaller than the prevailing wave¬ 
lengths, can produce significant oscillatory fluid flow, generally referred to as wave-induced fluid flow, 
between the different regions composing the heterogeneous medium (e.g., Muller et al. [2010 ). Due to 
the differing elastic compliances of the various regions, the stresses associated with the seismic per¬ 
turbation produce a pore fluid pressure gradient and, consequently, fluid flow. The energy dissipation 
related to this phenomenon is considered to be one of the most common and important seismic atten¬ 
uation mechanisms in the shallower parts of the crust (e.g., |Miiller et al.|2010| ). Given that the amount 
of fluid flow produced by this phenomenon can be significant, a potentially important interface-type 
seismoelectric signal is also expected to arise. However, the nature and importance of the seismo¬ 
electric effects related to mesoscopic heterogeneities remain largely unexplored. An interesting study 
of this kind was presented by Haartsen & Pride (19971 who modeled the seismoelectric response of 
a single sand layer having a thickness much smaller than the predominant seismic wavelength and 
embedded in a less permeable medium. They observed that while the seismic amplitudes recorded at 
the surface were very small due to destructive interferences, the converted EM amplitudes were sig¬ 
nificantly enhanced compared to the case of a thick sand layer. More recently,IJougnot et al.|(2013l 
proposed a numerical framework to study the seismoelectric response of a rock sample containing 
mesoscopic fractures subjected to an oscillatory compressibility test (Rubino et al7|[2009[ 20131 and 


found seismoelectric signals that would be measurable for typical laboratory setups (e.g., Tisato & 
Quintal! 20131. These findings are important not only as the seismoelectric responses of most geo¬ 
logical environments arc expected to contain a component related to fluid flow at mesoscopic scale, 
but also because they open the perspective of developing seismoelectric spectroscopy as a novel lab¬ 
oratory method for characterizing heterogeneous porous media. Jougnot et ah] (2013 ) suggest that a 
better understanding of the role played by mesoscopic heterogeneities could help to improve some of 
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the practical aspects of the seismoelectric method, such as the notoriously high noise levels generally 
observed in the measurements. 

In this work, we present an analytical approach to study the origin of the seismoelectric signal 
produced by mesoscopic heterogeneities. We first derive a general analytical solution for the seismo¬ 
electric response of a homogeneous rectangular rock sample containing a horizontal layer at its center 
and subjected to an oscillatory compressibility test. Following Brajanovski et al. (20051, we then adapt 
the analytical solution to the particular case of a sample containing a central layer having thickness and 
dry-frame elastic moduli tending to zero in conjunction with porosity tending to one. This particular 
solution represents the response of a sample permeated by a horizontal fracture at its center. Finally, 
we employ these two solutions to explore the roles played by different petrophysical and structural 
properties in the seismoelectric signatures of heterogeneous rocks. This analysis may, in turn, help to 
identify which parameters could be retrieved from seismoelectric measurements. 


2 METHODOLOGY 

2.1 Oscillatory compressibility test 


We consider a rectangular fluid-saturated porous rock sample containing a horizontal layer located 
at its center (Fig. |T]). The two regions embedding this mesoscopic heterogeneity are identical. For 
simplicity, we choose the center of the sample as the origin of the z- axis and, therefore, the positions 
of the upper and lower boundaries of the heterogeneity are z = L\ and z = —L\, respectively. The 
thicknesses of the two embedding regions constituting the host rock are L 2 and the total thickness of 
the sample is L = 2{L\ + L 2 ). 

The sample is subjected to a time-harmonic compression of the form /\Pe lujt at its top bound¬ 


ary, with 1 0 being the angular frequency, t the time, and i = (Rubino et al. 2009 2013). We 
impose that the solid phase is not allowed to move on the bottom boundary, nor to have horizontal dis¬ 
placements at the lateral boundaries and the pore fluid cannot flow into or out of the sample. For the 
considered geometry and boundary conditions, the problem to solve is ID. We consider frequencies ui 
smaller than Biot’s critical frequency lob (e.g., Bi ot|1962| ) 


0 , <fa 

lo b = 27 xJb =-, 

K f J f 


( 1 ) 


where <j) is the porosity, k the permeability, pj the density of the pore fluid, and // the fluid viscosity. In 


this frequency range, viscous boundary layer effects are negligible and thus we can solve Biot’s (1941) 


consolidation equations to obtain the response of the sample. In the ID space-frequency domain these 
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equations can be written as 


dr „ 

~dz ~ ' ’ 

t] dpt 

IU-W = — 

K OZ 


( 2 ) 


(3) 


where r and pf denote the total stress and the fluid pressure, respectively, whereas w is the relative 
fluid-solid displacement. Equation ([2]) represents the stress equilibrium within the sample, whereas 
Eq. ([3]) is Darcy’s law. These two equations are coupled through the ID stress-strain relations 


rr du „, dw 
T = H— + aM—, 

OZ OZ 


r du 


dw 


Pf = - aM - - M —. 


(4) 


(5) 


In these equations, u denotes the average displacement of the solid phase and the coefficients H, M 
and a arc given by 


M = 


H = A + 2/jl, 
a — <f> f 


K, 


K 


f J 


-l 


a = 1 — 


K 

lv m 

~Kl' 


( 6 ) 

(7) 

(B) 


where K m , K s , and Kf are the bulk moduli of the solid matrix, solid grains, and fluid phase, respec¬ 
tively. Moreover, fi is the shear modulus of the bulk material, which is equal to that of the dry matrix, 
and A is the Lame constant, given by 


A — Km T Klo? — —/x. 


For a homogeneous medium, combining Eqs. ([2]) and Q, as well as Eqs. (|3]) and (|5]) leads to 

d 2 u aM d 2 w 

dz 2 H dz 2 ’ 


(9) 


( 10 ) 


and 


V „ r 9 2 u , d 2 w 

mjd-w = aM—r + M ——, 
k dz z dz z 


respectively. Next, substituting Eq. ( fT0| ) into Eq. ( [IT] ), results in 

d 2 w 


iuw = D 


dz 2 


Equation ( [T2| ) is a diffusion equation, with the diffusivity D given by 

k N 


D = 


V 


( 11 ) 


( 12 ) 


(13) 


where N = M — a 2 M 2 /H. 
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Figure 1. Schematic representation of a rectangular porous rock sample containing a horizontal layer (white 
rectangle) embedded between two identical regions (gray upper and lower rectangles). The sample is subjected 
to an oscillatory compressibility test. 


The spatial scale at which wave-induced fluid flow is significant is determined by the diffusion 
length 

L d = yfDju. (14) 

From Eqs. and ( fl~4| > it is clear that this length increases with increasing permeability and de¬ 
creasing viscosity. When the diffusion length is of similar size as the characteristic size of the hetero¬ 
geneities, ameso, a characteristic frequency can be defined as 

D 


fc 


27r al 


(15) 


which represents the limit between two types of mechanical behaviors in response to the oscillatory 
compression. For frequencies / <C f c , the diffusion lengths are much larger than the typical size of 
the heterogeneities. Correspondingly, there will be enough time during each oscillatory half-cycle for 
the pore fluid pressure to equilibrate at a common value. Thus, this low-frequency regime represents 
a relaxed state. For frequencies / f c , the diffusion lengths are very small compared to the size of 
the heterogeneities and, hence, there is no time for interaction between the pore fluids of the different 
parts of the rock. In this case, the pore pressure is approximately constant within each heterogeneity 
and, consequently, the high-frequency regime is associated with an unrelaxed state. For intermediate 
frequencies, / ~ f c , significant fluid flow can be induced by an oscillatory stress field (e.g., Muller 


et al. |2010| >. From Eqs. (p~3|) and <fT3|) we notice that the frequency range where significant fluid flow 


can occur shifts towards lower frequencies for decreasing permeability, increasing fluid viscosity, or 
increasing size of the heterogeneities. 
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The general solution of Eq. ( [T2| ) for a homogeneous medium is given by 

w{z) = Ae~ kz + Be kz , 


(16) 


where 



(17) 


and A and B being complex-valued constants. From Eq. ([4]), given that r is spatially constant in virtue 
of Eq. Q, the solid displacement can be expressed as 


u(z) = —/3w(z) + ^—z + C, 


(18) 


where /3 = aM/H is the ID Skempton coefficient (Appendix A) and C is an additional complex¬ 
valued constant. Eqs. © to ( |18| ) then constitute the general solutions of Eqs. ([2]) to ([ 5 ]) for homoge¬ 
neous media. 

Since the generation of the seismoelectric signal is produced by the relative fluid velocity field, 
iujw, we seek an analytical expression for this parameter. The solid displacement u always appears 
as a derivate with respect to z in Eqs. to ([ 5 ]). This means that if (w(z), u(z )) constitute a solution 
for such equations, then ( w(z),u(z ) + x), with x being a constant, is also a solution. Hence, the 
solution for w is independent of the solid displacement value imposed at the bottom boundary of the 
sample. Correspondingly, without changing the solution for w, we can modify the boundary condition 
imposed on u at the bottom boundary to produce u = 0 at z = 0. Under this condition, and since r is 
constant within the sample (Eq. ([2])), it is clear that the geometry and the imposed boundary conditions 
are symmetrical with respect to z = 0 and, therefore, the resulting fluid flow profile should also be 
symmetrical with respect to the center of the sample. Hence, we can simply solve Eq. ( |T?| ) in the upper 
half of the sample shown in Fig. [I] under the following boundary conditions 


t = —A P, z = LI 2, (19) 

w = 0, z = L/2, (20) 

w = 0, z = 0, (21) 

u = 0, z = 0. (22) 


According to Eqs. ( jT6| ) and ( p~8] >, the general solutions in the upper half of the sample can be expressed 
as 


Wj{z) = A je ~ k i z + B je k * z , 

A P 

uj(z ) = —(3jWj(z) ——z + Cj, 

tl n 


where r has been replaced by the imposed stress 


(23) 

7 = 1,2, (24) 

—A P. The subscripts 1 and 2 denote the parameters 
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of the lower (0 < z < L\) and upper {L\ < z < L/2) layers composing the upper half of the sample, 
respectively. The six unknowns Aj, Bj, and C } (j = 1, 2) are determined by imposing the conditions 


given by Eqs. (20) to (22 1 and the continuity of u, w, and pj at the interface z = L\. Taking these six 
conditions into account and the fact that the relative fluid displacement w is an odd function, it can be 
shown that 


w(z) = < 


sgn{z)A\ (e fcl l 2 l — e fcl l z l) , 0 < \z\ < L\, 

sgn(z)A 2 [e.~ k2 ^ — , L\ < \z\ < L/ 2, 

where sgn is the sign function. The parameters A\ and An are given by 


(25) 


Ai = [e~ klLl -e klLl 


-l 


AP(/3i-/3 2 ) 


An = ' « _fc2Ll 


— e 


-fc 2 (L-Li) 


Njkj coth(kjLj) ’ 

-i AP(ft-fe) 

E?=i N i k j coth(kjLj) 


(26) 


(27) 


In the following subsection, the above expressions for the relative fluid displacement are used to 
infer the distribution of the electrical potential within the sample in response to the applied oscillatory 
compression. 


2.2 General solution: seismoelectric response of a rock sample containing a horizontal layer 

When a relative displacement between the pore fluid and the solid frame occurs in response to the 
applied oscillatory compression, a drag on the electrical excess charges of the EDL takes place. This, 
in turn, generates a source or streaming current density j s . Since the distributions of both the excess 
charge and the microscopic relative velocity of the pore fluid are highly dependent on their distance to 
the mineral grains, not all the excess charge is dragged at the same velocity by the flow. Correspond¬ 
ingly, an effective excess charge density Q^ smaller than the total excess charge density Q v has to be 
employed (e.g., Jougnot ct al. |2Q 1 2\ Revil & Mahard lka|2013 1. In the considered ID case, the source 
current density takes the form (e.g., Jardani et al.|2010t Jougnot et al.|2013 1 


js = Qfiuw, 


(28) 


where iurw is the relative fluid velocity. The effective excess charge formulation, which allows us to ex¬ 
plicitly express the role played by the relative fluid velocity in the source current density generation, is 
equivalent to the electrokinetic coupling coefficient formulation commonly used in the seismoelectric 
literature (e.g., |Pride|]l994 ). The relationship between the effective excess charge and the electroki¬ 
netic coupling coefficient can be found in many works such as Revil & Leroy ( 2004] ), Jougnot et al. 
|20l2| ), or [Revil & Mahardika| ( |2013| ). 
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In the absence of an external current density, the ID electrical potential p in response to a given 
source current density satisfies (Sill|19831 


d f dp 
dz \ dz 


dj s 

dz ’ 


(29) 


where o denotes the electrical conductivity, which strongly depends on the saturating fluid as well as 
on textural properties of the medium, such as porosity and tortuosity. 

The approach presented above is based on the fact that the electrical potential field has a negligible 
effect on the fluid flow pattern and, therefore, the seismic and electrical problems can be assumed to 
be decoupled. This assumption is made in most seismoelectric studies for materials similar to the ones 
considered in this work (e.g., Haines & Pride|2006] Guan & Hu|2008[ Zyserman et al.|20T0 i. 

In order to obtain the seismoelectric response of the sample shown in Fig. [T] we must solve Eqs. 
( [28] ) and ( [29} for the relative fluid displacement given by Eq. ([25]). The general solution for the elec¬ 
trical potential p is then given by 


p(z) = < 



^2 

<72 

/C2 

1 

At 

cr 1 

ki 

iuQf 2 

A2 


k 2 


(e k2Z + e~ k ^ L+z ^) +Rfz + Sf, 


i (e- k ' z + e k ' z ) +R 1 z + S 1 , 

(e~ k2Z + g-fofi-z)) + R+ z + s+, 


—L/2 < ^ < -Li, 
—L\ < z < L i, 
L 1 <z< L/2. 


(30) 


The parameters R 2 , S 2 , R\, S\, R 2 and S 2 are complex-valued constants that can be obtained by 
imposing the following boundary conditions 


= 0 
dz 

P = 0, 


2 = -L/2, L/2, 
z = L/2. 


(31) 

(32) 


Equation ( |Ti~[ ) states that the rock sample is electrically insulated at its bottom and top boundaries, 
whereas Eq. ([32]) indicates that the top boundary is the zero reference for the electrical potential. 


To obtain an additional condition, we integrate Eq. (29 1 in the upper half of the sample 

L/2 


a 


dp 

dz 


_ | L /2 
— 2s lo 


(33) 


Since w = 0 in both the top boundary and the center of the sample, the right-hand side of Eq. (33 1 is 
zero. Using Eq. (f3l[) then gives 


dp 

= 0 , 2 = 0 . 
dz 


(34) 


The four boundary conditions stated in Eqs. <[3T}, ©, and ( [34} , together with the continuity of p(z) 
at z = ±T], provide six relations that allow us to determine the complex-valued constants of Eq. ([30}. 
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V(z) = < 


By solving this linear system of equations, we obtain 

{e~ kl]z \ + e fel|z| ) + Si, 

-^2^ (e-HA + e -fc 2 (£-M)) + s 2 , 

where S\ and S 2 are given by 

s = iuj &Z m ( c ~k lL 


0 < \z\ < L\, 
Li < \z\ < L/2, 


01 h 


+ e 


kiLi 


i“Q%,2 A 2 / _ k2Li 
a 2 k 2 V 

2 ivQl G 2 A 2 k l / 2 

=-e 

cr 2 k 2 


+ e 


-k 2 {L-L i) 


- 2e 


( 35 ) 


-k 2 L/ 2^ 


(36) 


(37) 


Equation ( [33] ), together with Eqs. ( [36] ) and (37 1 , constitute the analytical solution describing the seis- 
moelectric response of a rock sample containing a central horizontal layer that is subjected to an 
oscillatory compressibility test (Fig. [Tj). 


2.3 Particular solution: seismoelectric response of a rock sample containing a horizontal 
fracture 

Fractures are present in virtually all geological formations, and they tend to control the overall hy¬ 
draulic and mechanical properties of these formations. This is why there is great interest in developing 
techniques to characterize fractured materials. Due to the very strong compressibility contrast typi¬ 
cally observed between fractures and the host rock, wave-induced fluid flow and, therefore, the corre¬ 
sponding seismoelectric effects, are expected to be significant in these environments. Indeed, Jougnot 
et al. ( (2013 l showed that under typical laboratory conditions the presence of fractures can produce 
measurable seismoelectric signals in response to the application of an oscillatory compressibility test. 
Here, we adapt the general analytical solution derived above to the particular case of a rock sample 
containing a horizontal fracture at its center. 


The poroelastic response of a fractured rock can be obtained in the framework of Biot’s (11962 1 
theory by representing the fractures as highly compliant and permeable heterogeneities embedded in 
a stiffer and less permeable host rock. In this sense, Brajanovski et al. (2005 [) employed an analytical 


solution of Biot’s (1962) equations for periodically varying coefficients to compute seismic attenu¬ 
ation and dispersion due to wave-induced fluid flow in fractured rocks. They considered the limit of 
very small values for the stiffnesses and apertures of the fractures in conjunction with very high poros¬ 
ity, which allowed for obtaining simple expressions for the corresponding effective complex-valued 
frequency-dependent plane wave modulus. 

Following the ideas of Brajanovski et al. ( [2005 l, we propose an analytical approach to obtain the 
seismoelectric response of a rock sample containing a horizontal fracture at its center. For this puipose, 
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the analytical solution obtained for a horizontal layer can be appropriately adapted by considering in 


Eq. (35 i that the aperture and porosity of the fracture satisfy h = 2L\ —y 0 and cff 1. Under these 
assumptions, the contribution of the fracture to wave-induced fluid flow can be significant only if the 
fracture stiffness also becomes very small. To take into account this inter-dependence, and following 
Brajanovski et ah](2005), we characterize the elastic properties of the drained fracture in terms of the 


shear compliance Zt and drained normal compliance Zy, which are given by 
Zt = lim — 

h ^° f 4 n 

h 

Zn = lim 


(38) 

(39) 


h ~*° Kh + \nln ’ 

where K/ n and f dn are the drained-frame bulk and shear moduli of the fracture, respectively. By taking 
the above limits, and using the fracture shear and normal compliances, Eq. ([35]) takes the form 


M ( e -kM + e -k 2 (L-\ Z I) _ 2e -k 2 (L/ 2)\ 

(T? ko \ ) 


cr 2 k 2 


(40) 


where 

A 2 = lim A 2 = - 2 -r- A^(l -fa) - (4 i) 

h ->o _ e -k 2 _|_ jv 2 fc 2 (1 + e _fc2L ) 

This analytical solution allows for computing the seismoelectric response of a rock sample containing 
a horizontal fracture at its center. Note that the only fracture parameter involved in these equations 
is the drained normal compliance Z ; .y and the only structural parameter is the total thickness of the 
sample L. 


3 RESULTS 

In this section, we employ the general and particular analytical solutions derived above to explore the 
roles played by different petrophysical and structural properties of heterogeneous porous rocks on the 
seismoelectric signals produced by oscillatory compressibility tests. In all cases, we assume that the 
pore fluid is water and that the solid grains consist of quartz (see Table [T]). 


3.1 Petrophysical relationships 

For the analysis based on the general solution, we consider a poroelastic model corresponding to a 
rectangular rock sample containing a horizontal layer at its center (Fig. [T]). Both the layer as well as 
the host rock correspond to clean sandstones, albeit with different porosities. To relate the porosity f 
to the permeability k, we use the Kozeny-Carman equation (e.g., Mavko et al. 2009) 


(i - <t>y 


ff 2 


k = b 


(42) 
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Table 1 . Material properties employed in the analysis considered in this work. 


Quartz grain bulk modulus, K s [GPa] 

Quartz grain shear modulus, fj, s [GPa] 

Water bulk modulus, Kf [GPa] 

Water viscosity, rj [Pa x s] 

Water electrical conductivity, a w [S m -1 ] 
Water density, p w \K.g m -3 ] 

37 

44 

2.25 

0.001 

0.01 

10 3 




Material 1 

Material 2 

Material 3 

Porosity, (j> 

0.05 

0.2 

0.4 

Dry rock bulk modulus, K m [GPa] 

31.47 

16.02 

2.88 

Dry rock shear modulus, p [GPa] 

37.42 

19.05 

3.42 

Permeability, n [mD] 

2.66 

240 

3410 

Electrical conductivity, a [S m _1 ] 

2.5 x 10" 5 

4 x 10~ 4 

1.6 x 10 -3 

Effective excess charge density, [C m -3 ] 

526.8 

13.13 

1.49 

Biot’s critical frequency, fg [Hz] 

2.99 x 10 6 

1.32 x 10 5 

1.8 x 10 4 


where b is a geometrical factor that depends on the tortuosity of the porous medium, and d the mean 
grain diameter. In this analysis, we take b = 0.003 and d = 8 x 10 -5 m. In addition to changes in 
permeability, porosity variations also imply changes in the mechanical properties. To link the porosity 
and the solid grain properties with the elastic moduli of the dry frame, we use the empirical model of 


Krief et al. (1990) 


K m = K s (1 - 

^ml^s 


d = 


K, 


(43) 

(44) 


where /r s is the shear modulus of the solid grains. 

The electrical conductivities cr of the considered clean sandstones are obtained using the empirical 


relationship proposed by Archie (1942) 


cr = </) m (T w , 


(45) 


where m = 2 is the cementation exponent and a w the electrical conductivity of the pore water. The 
remaining electrical parameter, is obtained by employing the empirical relationship proposed by 


Jardani et al. (2007) 


log (Qf) = -9.2349 - 0.8219 log(/e), (46) 

where k and are in units of nr 2 and C/m 3 , respectively. Below Biot’s critical frequency the ef- 
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(a) (b) 




Figure 2. Comparison between (a) the real and (b) the imaginary parts of the electrical potential of a rock 
sample containing a horizontal layer at its center based on the analytical solution derived in this study (solid 
lines) and the numerical framework by Jougnot et al. ( 2013| > (circles) as functions of position. The different 
curves correspond to three different frequencies of the applied harmonic compression. The imaginary part is 
related to the lag between the applied oscillatory compression and the resulting electrical potential distribution. 


fective excess charge density is similar to the one at zero frequency (e.g., Tardif et al. 20ll]iRevil & 
jMahardika||2013j ) and, hence, boundary layer effects can be neglected in the test cases considered in 
the following. Even if we consider simple petrophysical relationships to link o and to porosity, 


these properties can also be determined independently in laboratory experiments (e.g., Jouniaux & 
Pozz ifi995||Suski et al.|2006j ). 


For the analysis of the particular solution, we consider a rectangular homogeneous sandstone 
containing a horizontal fracture at its center. Equations (42 1 to ( |46| ) are employed to determine the 
petrophysical properties of the host rock, while the parameter Zn characterizes those of the fracture. 


3.2 General solution analysis 

3.2.1 Compliant, high-permeability layer 

We first consider a rock sample with a vertical sidelength of 20 cm composed of a stiff, low-permeability 
host rock with a porosity of 0.05 (material 1 in Table [T]), permeated at its center by a compliant, high- 
permeability horizontal layer with a thickness of 6 cm and a porosity of 0.4 (material 3 in Table[l]). The 
sample is subjected to a harmonic compression of amplitude A P = 1 kPa and frequencies ranging 
from 1 to 10 4 Hz. Before analyzing the seismoelectric response of this sample, we show in Fig. [5] a 
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Figure 3. (a) Relative fluid velocity, (b) electrical source current density, and (c) amplitude and (d) phase of 
the electrical potential corresponding to a rectangular, stiff, low-permeability host rock containing a compliant, 
high-permeability horizontal layer at its center. The porosity of the layer is 0.4 (material 3 in Table [TJ, whereas 
that of the host rock is 0.05 (material 1 in Table[I](. In all cases, the panels show the parameters along the ::-axis 
as functions of the frequency of the applied harmonic compression. For visualization purposes, we indicate the 
boundaries of the layer using white lines. 


comparison between the electrical potential calculated with the analytical solution described in this 
work and the solution obtained using the numerical framework proposed by Jougnot et al. (2013). 
Due to the symmetry of the solution, we only show the response in the upper half of the sample. We 
observe excellent agreement between the electrical potential curves obtained using the analytical and 
numerical methodologies. 

Returning to the analytical solution for the compliant layer. Fig. |3ja) shows the amplitude of the 
resulting relative fluid velocity along the 2 -axis of the sample as a function of frequency. Due to the 
strong contrast between the compressibilities of the two materials, significant relative fluid velocities 
arise in both the host rock and the layer. There is a frequency at which the spatial extension of host rock 
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affected by fluid flow is at its maximum, which we denote as the peak frequency. The characteristic 
frequency f c given by Eq. < ]T3| ) is indicative of this frequency, since for such a value the diffusion length 
in the host rock is comparable to the characteristic size of the region. For lower frequencies, fluid 
velocities tend to become negligible. For higher frequencies, the regions of the host rock affected by 
fluid flow tend to only include the immediate vicinity of the boundaries of the layer, but the magnitude 
of the relative fluid velocity is important and increases with frequency. Relative fluid velocity is also 
important inside the layer, mainly in the vicinity of the boundaries. Since the permeability of this 
material is much higher, the corresponding peak frequency occurs at higher frequencies, as suggested 
by Eq. £3) and (B). 


A significant current density (Fig. 0 b)) prevails in the host rock due to the relative fluid velocity 
field (Fig. [3ja)) produced by the compression. Moreover, the source current density clearly depends 
on the frequency, a relation that arises from the frequency-dependence of the induced fluid flow. The 
maximum current densities occur at the contacts between the two materials, where the relative fluid 
velocity is also maximum. Within the layer, even though significant fluid flow also takes place, the 
resulting source current density is very small, since the effective excess charge is much smaller in this 
material characterized by a much higher permeability (Table 0 Eq. (|46|)). 


Significant electrical potential differences (Fig. 0c)), well above the ~ 0.01 mV threshold of lab¬ 
oratory experiments (e.g., Zhu & Toksoz|2005| Schakel ct al.|201 2), arise in response to the oscillatory 


compression. These results are consistent with those of Jougnot et al. ( 2013) for fractured rocks and 
points to the importance of wave-induced fluid flow effects on seismoelectric signals in the presence 
of porosity variations. In the host rock, the frequency-dependence of the electrical potential distribu¬ 
tion is very similar to that of the current density. Moreover, the region in the host rock characterized 
by significant values of electrical potential has its maximum spatial extension at the corresponding 
peak frequency. Inside the layer the behaviour is different. For each frequency, the amplitude of the 
electrical potential is rather constant. This happens because in this high-permeability material the ef¬ 
fective excess charge is small while the electrical conductivity is large, which implies that the electrical 
potential is reduced to a near-constant value, ,S'| (Eq. (|35|)). Because the electrical potential is contin¬ 
uous, this constant corresponds to the value of the electrical potential at the contact between the two 
materials. 


The resulting electrical potential is not only characterized by its amplitude but also by its phase. 
For relatively low frequencies, the phase 9 of the electrical potential remains constant, while for high 
frequencies it shows rapid changes within the host rock (Fig. [3jd)). Inside the layer 9 remains constant, 
in agreement with the behaviour observed for the amplitude of the electrical potential in this region 
(Fig. 0c)). 
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Figure 4. (a) Relative fluid velocity, (b) electrical source current density, and (c) amplitude and (d) phase of 
the electrical potential corresponding to a rectangular, compliant, high-permeability host rock containing a stiff, 
low-permeability horizontal layer at its center. The porosity of the layer is 0.05 (material 1 of Table[l]i, whereas 
that of the host rock is 0.4 (material 3 in Table[l]l. In all cases, the panels show the parameters along the ::-axis 
and as functions of the frequency of the applied harmonic compression. For visualization purposes, we indicate 
the boundaries of the layer using white lines. 


3.2.2 Stiff, low-permeability layer 

We now repeat the preceding analysis, but with the properties of the host rock and of the layer inter¬ 
changed. That is, we assume that the layer and the host rock have porosities of 0.05 (material 1) and 
0.4 (material 3), respectively, with the corresponding petrophysical properties given in Table [T] Due 
to the strong compressibility contrast between the two materials, significant fluid flow arises in the 
vicinity of the contact zones between the host rock and the layer (Fig. |4ja)). Moreover, we observe 
that the peak frequency corresponding to the stiff, low-permeability material is higher compared with 
the previous situation (Fig. |3ja)). The reason for this is that the thickness of the stiff material is now 
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Figure 5. Electrical potential amplitude \ip\ for a layer having a thickness of 2L-\ = 12 cm, that is, twice the 
value considered in the cases shown in Figs. [3] and [4] Panel (a) corresponds to a compliant, high-permeability 
layer (material 3 in Table [TJ embedded in a stiff, low-permeability host rock (material 1 in Table [TJ, whereas 
panel (b) shows the results obtained when the two materials are interchanged. For visualization purposes, we 
indicate the boundaries of the layer using white lines. 

smaller and, therefore, as dictated by Eq. & the characteristic frequency shifts towards higher fre¬ 
quencies compared to the situation depicted in Fig. j3ja). We also observe that for lower frequencies 
the fluid velocities tend to be negligible. For higher frequencies, significant fluid velocities prevail in 
smaller spatial extensions of the layer but their magnitudes increase with frequency. In the compli¬ 
ant host rock, the peak frequency is higher because this material is much more permeable (Eq. ©). 
which explains the differing behaviour of the fluid velocity as compared with that observed in the stiff 
material. 

The highest magnitudes of the electrical source current density are now concentrated inside the 
low permeability layer (Fig.|4jb)), which is characterized by high effective excess charges. As a result, 
and due to the imposed boundary conditions, when the layer is much stiffer and less permeable than 
the host rock, the electrical potential has a significant amplitude only inside the layer (Fig. |4jc)). 
Moreover, the electrical potential amplitude is frequency dependent, with a maximum at the centre of 
the layer and for a frequency close to the peak frequency (Fig. [4ja)). With respect to the phase of the 
electrical potential, significant phase changes within the sample take place at very high frequencies 
and in the host rock (Fig. [4jd)). 

3.2.3 Sensitivity to the thickness of the layer 

To analyze the role played by the thickness of the layer, we now consider 2Li = 12 cm for this region, 
that is, twice the original value, but leave the overall model setup and total thickness of the sample 
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otherwise unchanged. The frequency- and space-dependence as well as the magnitude of the electrical 
potential for the compliant and stiff layer models, shown in Figs. |5ja) and [5jb), are similar to the 
ones shown in Figs. |3jc) and |4jc), respectively. Flowever, the maxima of the electrical potential do 
not occur at the same frequencies. This shift in frequency can be explained by the dependence of 
the characteristic frequency f c , which indicates the frequency at which maximum penetration of fluid 


flow takes place, on the characteristic size of the involved medium. Indeed, Eq. (151 dictates that f c is 
inversely proportional to the characteristic size of the region where fluid flow occurs. When the layer 
is more compliant and permeable than the host rock, the electrical potential is mainly produced by 
fluid flow in the host rock. As the thickness of the layer increases with respect to the situation depicted 
in Fig. [3jc), the characteristic size of the host rock is reduced. This reduction, in turn, produces an 
increase of the characteristic frequency f c , which can be verified by comparing the frequency ranges 
where maximum electrical potential occurs in Fig. [3|c) and Fig. [5Ja). Conversely, when the layer is 
stiffer and less permeable than the host rock, the electrical potential is produced by fluid flow inside 
the layer. Therefore, as the thickness of this region increases in the situation shown in Fig. [5jb), the 
characteristic frequency f c is reduced, in agreement with Eq. S Correspondingly, the frequency 
range where the electrical potential is maximum shifts towards lower frequencies, as can be verified 
by comparing Fig. |4|c) and Fig. |5jb). This result therefore indicates that the seismoelectric signal is 
sensitive to the layer thickness, which largely controls the frequency range where the maximum signal 


occurs. 


3.2.4 Sensitivity to the compressibility contrast between layer and host rock 

Since the amount of fluid flow scales with the compressibility contrast between background and het¬ 
erogeneities (e.g., |Miiller et al.|2010 i, the porosity contrast between the layer and the host rock is ex¬ 
pected to have a major influence in the resulting seismoelectric response. To analyze this, we present 
in Fig. [6] the amplitude of the electrical potential as a function of frequency and space for the same 
geometries used in the cases shown in Fig. [5]but considering a reduced porosity contrast. We now use 
porosities of 0.05 (material 1 in Table [TJ and 0.2 (material 2 in Table |T[) for the stiff and compliant 
materials, respectively. By comparing Figs. [5] and Fig. [6] we can verify that the considered change 
in porosity contrast does not significantly affect the general frequency-dependence and spatial distri¬ 
bution of the electrical potential. This is because the geometry as well as the properties of the stiff 
material, which is the region where most of the current density arises, are similar. Conversely, there is 
an important decrease of the amplitude of the electrical potential. This effect is expected, as a reduction 
of porosity contrast implies a decrease of compressibility contrast. The latter causes a reduction of the 
wave-induced fluid flow (e.g., Muller et al.|2010i and, thus, of the source current density. These results 
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Figure 6. Electrical potential amplitude \ip\ for a rock sample containing a horizontal layer of thickness 2L \ = 
12 cm and characterized by a reduced porosity contrast as compared with the situation depicted in Fig. [5] Panel 
(a) corresponds to porosity values of 0.2 (material 2 in Table[l]i and 0.05 (material 1 in Table[I]( for the layer and 
host rock, respectively, whereas panel (b) shows the results obtained when the two materials are interchanged. 
For visualization purposes, we indicate the boundaries of the layer using white lines. 


illustrate that the magnitude of the electrical potential may contain information on the compressibility 
contrast between the layer and the host rock. 


3.2.5 Asymptotic analysis 


The previous parametric analysis suggests that key mechanical, structural and hydraulic properties 
of the central layer and of the host rock affect the seismoelectric response. In particular', while the 
compressibility contrast between the materials controls the magnitude of the resulting seismoelectric 
signal, the thickness and the permeability of the less permeable region determine the frequency range 
at which the response takes its maximum value. 

In order to provide a formal description of this dependence and to further explore the analytical 
expression given by Eq. ( |35j ), we develop the low- and high-frequency asymptotic behaviors of the 
solution. Since from the previous analysis we conclude that, for a given frequency, the electrical po¬ 
tential is at its maximum at z = 0, we derive these limits at this particular place of the sample. In 
Appendix B, we show that the low-frequency asymptote of the electrical potential at the center of the 
sample is given by 
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whereas the high-frequency asymptote is 




V f/o(u,) = -AP(A-ft) 


V 2 

^j=t V k j 


(48) 


The latter expression is the high-frequency limit of the solution obtained based on the premise that fre¬ 


quency is lower than Biot’s critical frequency. In other words, Eq. (|48]) represents the high-frequency 
limit of the low-frequency-domain solution. 

Figure [ 7 ] shows the amplitude of the electrical potential at the center of the sample as a function 
of frequency, for both the compliant- and stiff-layer cases considered in Figs. [3] and [4] together with 
the corresponding low- and high-frequency asymptotes. There is very good agreement between the 
general solutions and the low- and high-frequency asymptotes. Moreover, this analysis shows the 
importance of the analytical solutions obtained in this work and, in particular', of these asymptotes, 
since they allow us to explicitly observe how the various petrophysical parameters contribute to the 
seismoelectric response. The high-frequency asymptote, which is constant, can be considered as an 
indicator of the maximum magnitude of the electrical potential. Therefore, Eq. ( j48] ) allows for a direct 
appreciation of how the different petrophysical properties control the amplitude of the seismoelectric 
response. This maximum amplitude is directly proportional to the applied stress, to a weighted average 
of the effective excess charges of the two media, and to the Skempton coefficient contrast between the 
host rock and the layer. The Skempton coefficient is the ratio of the fluid pressure increase to the 
corresponding applied stress for undrained conditions (e.g.,|Wang|2000i. Correspondingly, the term 


APfa-fo), 


(49) 


appearing in Eq. ( [48] ) constitutes a measure of the fluid pressure gradient induced across the inter¬ 
faces separating the host rock and the layer. It is thus reasonable that the maximum amplitude of the 
seismoelectric response turns out to be directly proportional to this parameter, since the fluid pressure 
gradient induced by the applied compression triggers the fluid flow and the corresponding seismoelec¬ 
tric signal. 

Employing Eq. ( |46| ), the weighted average of the effective excess charges included in Eq. ( [48] ) 
satisfies 
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If the permeability contrast between the two media is strong, the numerator and denominator of the 
right-hand side of Eq. ([50]) are controlled by the term corresponding to the material of lower perme- 
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Figure 7. Amplitude of the electrical potential in the center of the sample (z = 0) as a function of frequency. 
Panel (a) corresponds to porosity values of 0.4 (material 3 in Table [TJ and 0.05 (material 1 in Table [7} for 
the layer and host rock, respectively, whereas panel (b) shows the results obtained when the two materials 
are interchanged. Red lines correspond to the exact solution (Eq. (f35]>), the blue ones to their low-frequency 
asymptotic solution (Eq. ([47])) and the pink ones to their high-frequency asymptotic solution (Eq. (|48|)). 


ability, that is, 
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where the subscript min denotes the material with the lowest permeability. Substituting this relation 
into Eq. ( [48] ), we obtain 


<P HF (uj) (X-AP (01 - 02 ) 


(52) 


Equation ([52]) suggests that the magnitude of the resulting electrical potential, in the case of strong 
permeability contrasts, is directly proportional to k^* 8 . This result also indicates that the hydraulic 
properties of the most permeable material and the thicknesses of the layer and of the host rock do not 
affect the magnitude of the seismoelectric response. This is in agreement with the results presented 
in this work, where similar values for the maximum amplitudes of the electrical potential were found 
when the materials of the host rock and the layer were interchanged. 

The frequency at which the low- and high-frequency asymptotes intercepts, f m t, can be used as 
an indicator of the frequency at which the electrical potential reaches its peak. This intersection can 
differ with respect to the true peak frequency by almost one order-of-magnitude (Fig. [7]). However, its 
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dependence on the mechanical, hydraulic, and structural properties of the heterogeneous sample are 
expected to also hold in the case of the true peak frequency. By requiring equality of Eqs. ( |47] ) and 
([48]). fi n t can be defined as 


fint 
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Hence, the frequency of the peak of the electrical potential does not depend on the amplitude of the 
applied pressure, A P. Interestingly, it does not depend on the Skempton coefficient contrast of the 
materials either. 

Considering the case of strong permeability contrasts and using Eq. ( [52] ) instead of Eq. ( |48j ) to 
obtain f int , yields 

j, Kmin Nj 

fint OC — - Z^T ■ ( 54 ) 
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3 =1 

This equation shows that the frequency of the peak of the electrical potential is directly proportional 
to Kmin, whereas the hydraulic properties of the most permeable regions do not affect this parameter. 
In addition, the structural characteristics of the sample do affect the location of the peak, which is 
inversely proportional to the thickness of the less permeable material. Also, the mechanical properties, 
scaled by the thicknesses of the corresponding layers, affect the peak frequency. 


3.3 Analysis of fractured media 


Fractures are usually composed of very compliant material and the resulting strong compressibility 
contrast with respect to the host rock is expected to favor particularly strong seismoelectric signals in 
response to the application of oscillatory compressions (Jougnot et al .[2013 1. In order to explore this in 
more detail, we estimate the seismoelectric responses of fractured media by making use of the model 
developed in Section 2.3. According to this model, in which the aperture of the fracture is considered 
negligible, the seismoelectric response is controlled by the properties of the host rock, the size of the 
sample, and the compliance of the fracture. The latter is represented by the drained normal compliance 
Z N (Eq. @). 

To study the behaviour of the seismoelectric response for different characteristics of fracture and 
host rock materials, we consider a rectangular rock sample with vertical sidelength of 20 cm containing 
a horizontal fracture at its center. As in the previous analysis, the sample is subjected to a harmonic 
compression with an amplitude A P = 1 kPa and frequencies ranging from 1 to 10 1 Hz. For the host 
rock, we consider a sandstone with a porosity of 0.05 (material 1 in Table [TJ, while, as proposed by 
Nakagawa & Schoenberg (2007 ), we use a drained normal compliance Zjy = 10 -8 m kPa 1 for the 
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Figure 8. Absolute value of the electrical potential \tp\ along the z-axis as a function of frequency for samples 
containing a horizontal fracture at their centers. Panel (a) corresponds to a host rock with a porosity o of 0.05 
(material 1 in Table [l} and a fracture drained normal compliance Z\r of 10 _8 m l<Pa 1 ; in panel (b) 0 = 0.4 
(material 3 in Table[l]) and Zpj = 10 _8 m kPa -1 ; in panel (c) <fi = 0.05 and Zm = 10 _6 m kPa^ 1 , and in panel 
(d) (f> = 0.4 and Zn = 10~ 6 m kPa” 1 . 


fracture. Very strong electrical potential differences, that could be measured using standard laboratory 
techniques, arise in response to the applied oscillatory compression (Fig. [8ja)). Moreover, the overall 
frequency- and space-dependence characteristics of the electrical potential are essentially the same as 
in the case of a rock sample containing a horizontal compliant and permeable layer. Indeed, there is 
also a frequency range at which both the magnitude and the spatial extension of the electrical potential 
are at a maximum. 

To study the role played by the characteristics of the host rock, we change its properties to those 
of a sandstone with a porosity of 0.4 (material 3 in Table [TJ. The corresponding increases of the 
permeability and of the compressibility of the host rock produce significant changes in the resulting 
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electrical potential (Fig. [8jb)). Indeed, the magnitude of the signal is significantly smaller than that 
shown in Fig.[8ja), which was evaluated using smaller values of permeability and compressibility of 
the host rock. This, in turn, indicates that, according to Eq. ( [52] ), the effect related to the reduction 
of compressibility contrast between fracture and host rock overcomes that related to the increase of 


the permeability of the host rock. In addition, as suggested by Eq. ( |54| ), the increase of the host rock 
permeability shifts the frequency range where the maximum value and maximum spatial penetration 
of the electrical signal occur towards higher values. 

The compliance of the fracture, represented by Zn, controls the compressibility contrast with 
regard to the host rock. This parameter is thus expected to play a major role in the seismoelectric 
response of fractured materials. To verify this, we show in Fig. [8jc) the electrical potential for a 
fractured sample having a host rock porosity of 0.05, that is, the same value employed to calculate 
the electrical potential represented in Fig. [8|a), but with a drained normal compliance Zy = 10 -6 m 
kPa -1 . The comparison between Figs, [sj a) and [sjcj shows that the overall frequency-dependence and 
spatial distribution of the electrical potential are very similar. However, an increase of fracture drained 
normal compliance implies an increase of compressibility contrast between fracture and host rock. 
This results in a significantly larger magnitude of the electrical potential. 

Finally, Fig.[8jd) shows the amplitude of the electrical potential for a host rock corresponding to a 
sandstone with a porosity of 0.4 and a fracture with a drained normal compliance Zjv = 10 _6 m kPa '. 
By comparing this electrical potential with the one represented in Fig. [8]», we see a significant reduc¬ 
tion of the magnitude of the signal together with a shift towards higher values of the peak frequency. 
The reduction in magnitude indicates that the increase of host rock compressibility related to the poros¬ 
ity change affects the compressibility contrast more significantly than the increase of compliance of 
the fracture. In addition, as dictated by Eq. the increase of host rock permeability related to the 
porosity variation shifts towards higher values the peak frequency of the electrical potential. 


4 DISCUSSION 


In accordance with Jougnot et al.| (20131, our results indicate that seismoelectric signals are highly 
sensitive to the presence of mesoscopic heterogeneities. The analytical solutions obtained in this work 
provide the explicit dependence of these signals on the structural, mechanical and hydraulic properties 
of the materials involved. In the following, we discuss some of the implications of this dependence. 

Assuming a laboratory experiment where electrodes are placed vertically along the sample shown 
in Fig. [T] our results suggest that one could identify the presence of a compliant and permeable layer 
as the region where the amplitude of the electrical potential remains constant, and of a stiff and less 
permeable layer as the region exhibiting a significant electrical potential gradient. Alternatively, the 
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center of the layer can be identified by finding the position of the maximum amplitude of the potential. 
In the case of fractured rocks, the drained normal compliance is a parameter of key importance 


for simulating seismic wave propagation in the framework of the linear slip theory (e.g., Schoenberg 


1980), but it is often poorly constrained. Our analytical solution for the seismoelectric response of a 


fractured rock could conceivably be employed to estimate the drained normal compliance of a fracture 
based on electrical potential measurements. First, the position of the fracture could be identified as the 
region where the maximum amplitude is measured and, then, the drained normal compliance could 


be estimated by inserting the measured value of electrical potential in Eq. (40). Additionally, the 
permeability of the host rock could also be retrieved, provided that measurements are made at several 
distances from the fracture or at different frequencies. 

The present analysis considers a rock sample subjected to an oscillatory compressibility test, 
which allows us to explore seismoelectric signals related to the generation of Biot’s slow waves at 
interfaces separating layers with differing mechanical properties. These results have interesting im¬ 
plications related to seismic wave propagation. The seismoelectric signals are expected to decay as 
the amplitude of the seismic wave progressively diminishes due to geometrical divergence, scattering, 
and attenuation effects. Therefore, in the case of a compressional seismic wave propagating through a 
periodic ID medium composed of an assembly of the sample shown in Fig. |T] similar to the situation 


studied by White et al. (1975) in the context of seismic attenuation due to wave-induced fluid flow, the 
electrical-potential contribution is expected to be similar to that shown in this work, but scaled accord¬ 
ing to the local stress applied by the passing seismic wave. A complete extension of the methodology 
presented in this work to seismic wave propagation is, however, not straightforward. In addition to 
scattering and amplitude distortion effects, further contributions, such as the effects of wavelength- 
scale relative fluid flow produced by the seismic wave, should be taken into account. The analysis of 
numerical simulations of seismic wave propagation through porous media containing mesoscopic het¬ 
erogeneities and of the resulting seismoelectric responses will be the subject of forthcoming studies. 

Among the implications of the explicit dependence of the electrical potential on rock parameters 
presented in this work, it is interesting to highlight that the magnitude of the resulting signal is directly 
proportional to the Skempton coefficient contrast between the involved materials, to the amplitude of 
the applied compression and, in the case of strong permeability contrasts, it is also highly influenced by 
the permeability of the less permeable material. In such a situation the peak frequency of the electrical 
potential is largely governed by the permeability and the thickness of the region containing the less 
permeable material. Any attempt to model the coherent noise observed in seismoelectric signals arising 
potentially from mesoscopic heterogeneities, as suggested by Jougnot et al. (2013), should care fully 
constrain the values of these parameters. Conversely, in the case of strong permeability contrasts the 
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hydraulic and electrical properties of the most permeable material do not need to be determined with 
high precision, as the seismoelectric signals turned out to be virtually insensitive to these parameters. 


5 CONCLUSIONS 

We have derived an analytical solution that describes the seismoelectric response of a rectangular rock 
sample containing a horizontal layer at its center that is subjected to an oscillatory compressibility 
experiment. The resulting solution was also adapted to compute the seismoelectric response of a rock 
sample containing a horizontal fracture at its center. 

The seismoelectric responses predicted by the analytical solutions were explored and allowed us 
to shed some light into their dependence on mechanical, hydraulic, and structural properties of the 
host rock and the layer or fracture. As a general result, we found that the seismoelectric signals pro¬ 
duced by oscillatory compressibility tests could be measured in standard laboratory experiments. The 
significant amplitudes of the seismoelectric signals point to the importance of considering mesoscopic 
effects when employing the seismoelectic method. Moreover, the analysis indicated that the maximum 
amplitude of the electrical potential is directly proportional to the applied stress, to the Skempton coef¬ 
ficient contrast between the host rock and the heterogeneity, and to a weighted average of the effective 
excess charge of the two materials. In presence of strong permeability variations we found that this 
weighted average is mainly controlled by the permeability of the less permeable material and that 
the frequency at which the maximum electrical potential prevails is governed by the permeability and 
thickness of the region containing such material. 

Although the analytical solutions considered in this work arc based on simple ID models, they 
constitute useful tools for exploring the details of the physical processes involved in the seismoelec¬ 
tric method in presence of mesoscopic heterogeneities. Moreover, they may open the possibility of 
retrieving key rock properties, such as permeability and fracture drained normal compliance, from 
seismoelectric measurements. 
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APPENDIX A: ID SKEMPTON COEFFICIENT 


The Skempton coefficient is defined as the ratio between the induced fluid pressure increase to the 
applied stress for undrained conditions (e.g., Wan g|2000 1 




C=o 


(A.l) 


where £ is the increment of fluid content, which in ID is given by C, = —dw/dz. This parameter 
measures how the applied stress is distributed between the solid frame and the pore fluid (Wang ]2000 ). 
Imposing the undrained condition £ = 0 in Eqs. (|4j) and ([5]), we get 
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The ID Skempton coefficient is thus obtained by taking the ratio between Eqs. ( |A.2[ ) and ( |A.3 1 
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APPENDIX B: ASYMPTOTIC ANALYSIS OF THE ELECTRICAL POTENTIAL AT THE 
CENTER OF THE SAMPLE 


Equation ( [35] ), together with the Eqs. ([36]) and ( |37[ ), constitute an analytical model for describing 
the seismoelectric response of a rock sample containing a horizontal layer that is subjected to an 
oscillatory compressibility test. Now we focus on the values of ip(z) at z = 0, that is, at the center of 
the sample. From Eq. (|35[), we get for the electrical potential 






fcl 


(B.l) 


Replacing the corresponding expressions for A\ and S\ (Eqs. (26) and (36) respectively) in (B.l) we 
obtain 

1 —cosh(kjLj) 


<p( o) = 


ipAPVh-MZUfk 


sinh(kjLj) 


£j=i N 3 k 3 coth(kjLj) 


(B.2) 


This equation predicts that the amplitude of the electrical potential is proportional to the Skempton 
coefficient contrast of the two involved materials and to the amplitude of the applied compression. 
However, the dependence of <p(0) on Q®®, oj and Lj and k j is not straightforward because of the 
hyperbolic functions involved. Developing the low- and high-frequency asymptotes of this equation 
may help understand the roles played by the different mechanical, hydraulic and structural properties 
on the seismoelectric response. 
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B1 Low-frequency asymptote 

Using first-order Taylor expansions the following approximations can be obtained 


1 — cosh(kjLj) 
sinh(kjLj ) 

kjLj coth(kjLj , 


- k i L v 


1. 


(B-3) 

(B.4) 


Eqs. (B.3 i and (B.4) are good approximations when kjLj —> 0, which occurs when cj —> 0. Replacing 


these first-order approximations in Eq. (B.2i, we get the low-frequency asymptote for the electrical 
potential in the center of the sample 


<&{<*) = 




£=i 


(B.5) 


B2 High-frequency asymptote 

In order to explore the high-frequency behaviour of </?(0), we use the following approximations which 
are valid when kjLj —t oo 


1 — cosh(kjLj) 
sinh(kjLj) 

coth(kjLj ) ~ 1. 


(B-6) 

(B.7) 


Replacing these expressions in (B.2 1 we obtain for the high-frequency asymptote 


-iuAP (/3i - &) £?_ i 

,HF / s _ (Tjkj 

t L z =0\U) — 


£j=i Njkj 

The expression for kj as a function of u can be obtained by combining Eqs.([T3|). (|T4|) and ([TT]) 
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(B.9) 


Then, replacing (B.9) in (B.8) one obtains 


L>z =oM = 


0 r)eff 

-a 

/ Nj 

Alj = 1 V K, 


(B. 10) 
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